Double heavy-quarkonium production from electron-positron annihilation in the 

Bethe-Salpeter formalism 
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In this work we evaluate the cross section of the process e + e~ — » J/ip n c at energy i/s w 10.6 GeV 
in the Bethe-Salpeter formalism. To simplify our calculation, the heavy quark limit is employed. 
Without taking the beyond-leading-order contribution(s) into account, the cross section calculated in 
this scenario is comparable with the experimental data. We also present our prediction for the cross 
section of double bottomonium production e + e~ — > T(lS)r]b for the energy range of yfs « (25-30) 
GeV which may be experimentally tested, even though there is no facility of this range available at 
present yet. 
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I. INTRODUCTION 

It is well known that there is a significant discrepancy 
between the experimental measurements [l], Q and the 
NRQCD predictions Hi] 

for the process e+e — > J/ip rj c 
at energy y/s « 10.6 GeV. To reduce the discrepancy 
many efforts have been made. For example, as discussed 
in Ref. Q, corrections from pure electromagnetic in- 
teractions are introduced into the non-relativistic QCD 
(NRQCD) factorization formalism; the next-to-leading- 
order contributions of strong interaction are taken into 
account in Refs. [6j, |7j; in Ref. [8|, the authors take into 
account corrections to the J/ leptonic decay width and 
the scale dependence of the leading-order prediction; etc. 

Indeed, NRQCD should work well while dealing with 
the processes where heavy flavors are involved, especially 
in this concerned reaction where only heavy flavors ex- 
ist. In the scenario of NRQCD an expansion is carried 
out with respect to powers of the relative three-velocity 
v which reflects the relativistic effects. Even though, in 
this reaction, the relativistic effects obviously are smaller 
than in the processes dominated by light flavors, they 
still may cause sizable contributions. Later works con- 
firm this viewpoint. One needs, therefore, to incorporate 
the relativistic corrections into the theoretical calcula- 
tions in an appropriate way. There are various ways to 
take into account the relativistic corrections. Ebert et al. 
carefully studied the corrections from relativistic effects 
and indicated that with the relativistic effects being 
properly handled, a prediction for the cross section which 
is consistent with the present experimental measurement 
can be obtained. 

Moreover, the factorization scheme should work in this 
case, because the production of J/tp + r\ c from e + e~~ via 
electromagnetic interaction can be regarded as a reversed 
Drell-Yan process where the factorization is proved. Thus 
one can factorize the hard and soft processes and then a 
convolution integration over the two parts results in the 
final amplitude. Therefore, one only needs to consider 



the relativistic corrections appearing in the soft part, i.e. 
in the hadronization process because the hard part is 
carried out in the framework of quantum field theory 
which is completely in the relativistic covariant form. 

To understand the experimental results, many authors 
have proposed various projects to improve the theoreti- 
cal framework in addition to NRQCD. For example, this 
problem was discussed with the method of perturbative 
QCD (pQCD) . In Ref. PH, the corrections of higher- 
twist wave functions were included in pQCD and the 
light-front quark model. 

Discussions given in Refs. Jf| @, HH, 03 suggest that to 
reduce the large discrepancy between the experimental 
results and the theoretical predictions based on NRQCD 
for the exclusive process with the final state which is 
composed of two charmonia, large next-to-leading-order 
(NLO) corrections may appear (the 'NLO' contribution 
is about 1.8 to 2.1 times of the leading-order one). In- 
cluding this large NLO contributions, their results are 
close to the lower bound set by the Babar and Belle col- 
laborations for the double-charmonia production. The 
authors also indicate that including the relativistic cor- 
rections can further enhance the estimated value. 

We would rather incline to believe that because of the 
validity of the factorization, perturbative calculation is 
suitable and, therefore, in our calculation we will ignore 
the contribution of the next-to-leading order corrections. 
Then one should expect that a larger correction may 
come from the soft part, i.e. the relativistic effect may 
be significant. 

As mentioned before, Ebert et al. considered such ef- 
fects . Alternatively, in this work, for properly incorpo- 
rating relativistic effects, we try to evaluate this exclu- 
sive process in the Bethe-Salpeter (BS) formalism [l3j . 
The BS equation is in principle established in the frame- 
work of relativistic quantum field theory, therefore, it is 
supposed to include all relativistic effects. Of course, to 
solve this equation in practice, one needs to adopt some 
approximations such as the instantaneous approximation 
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where part of the relativistic effects are lost. However, 
in many cases, such loss is not serious. Indeed, the BS 
formalism is suitable for studying the binding systems 
composed of two heavy charm and anti-charm quarks (or 
bottom flavors). The transition amplitudes can be ob- 
tained, in a natural way, to be an overlap integration of 
BS wave functions (see e.g. Ref. (l4|). In order to sim- 
plify the calculation, we will further impose the heavy 
quark limit [HI, Oil [13, i-e. all the 1/rriQ corrections are 
neglected in the calculation in this paper. Under this 
limit, our result shows that, at the leading order of a s , 
the theoretical prediction is comparable with the exper- 
imental result P, 0. We attribute this mainly to the 
inclusion of relativistic effects in our formalism. We will 
come back to give more discussion on this point in the 
last section. 

The remainder of this paper is organized as follows. In 
Sec. mi we will study the BS equations for the vector and 
pseudo-scalar quarkoniums. In Sec. IIII1 we will calculate 



the cross section of the process e H 



J/tjj rj c in the BS 



formalism. In this section, we will also give the prediction 
for the much smaller cross section of the exclusive process 



with double bottomonium production, e" 1 



T(lS)»fe. 



Sec. HVlis reserved for our conclusions and discussions. 



II. BS EQUATIONS FOR HEAVY 
QUARKONIUMS 



following way [2C 



1 



(1 + 0/2 



-niQ + ie p e + Eo/2 + ie 
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(3) 
(4) 



where we have defined the binding energy E = M— 2toq, 
v = P/M is the 'four- velocity' of the meson, M is the 
mass of the meson andp f = v-p. With this simplification, 
the BS equation @ becomes 



x P (p) 



(1 + 0/2 
+ E /2 + ie 

(i - m 

p, - E /2 - ie 



^V P {p,k) Xp {k) 



(5) 



From this expression, one can easily see that the BS wave 
function \ P satisfies = X P an d X P $ — — X P - Then, 
in the heavy quark limit, similar to the case for the di- 
quark system studied in Ref. [2l|, we have the following 
very simple parameterizations for the BS wave functions 
of vector and pseudoscalar QQ mesons (which are de- 
noted by the subscripts " a" and " b" , respectively) : 



Xp s) (p) = (l + ^)/ s) M a f a (p), J 



x P » = (i + jOt&aWp), j 



PC 



PC 



= 1- 
o-+. 



(6) 
(7) 



The BS wave function for a meson which is composed 
of a quark and an anti-quark is defined as 



X P (xi,X 2 ) a /3 



V3 



(1) 



where P is the momentum of the meson, ip* a (xi) and 

ipg{x2) are the field operators of the quark and anti- 
quark, respectively, a, (3 are Dirac spinor indices, and 
i,j denote the color indices. In momentum space the 
BS equation for the wave function x P { x ii ^2)0/3 can be 
written as (see Refs. [TH, [l9| for example) 



x P (p) 



mi + *e 
i 

• — TO 2 + 



(2) 



where p, k represent the relative momenta between the 
quark and anti-quark, nii(pi) and m 2 {p2) are the masses 
(momenta) of the quark and anti-quark, respectively, and 
Vp (p, k) is the kernel. The spinor indices are suppressed 
for simplicity. For the heavy quarkonium QQ (Q = c, b) 
studied in this paper, we have ui\ — m 2 = tuq and then 
p x = P/2 + p and p 2 = -P/2 + p. 

To simplify the calculation we will take the heavy 
quark limit throughout this paper. In this limit, the 
propagators of the heavy quark and heavy anti-quark 
(S(pi) and S(p 2 ) 7 respectively) can be simplified in the 



where e^> is the polarization vector of the vector quarko- 
nium which is orthogonal to the velocity, • v = 0, f a 
and ff, are scalar functions of the relative momentum p. 

In this work, following the standard procedure for solv- 
ing the BS equation we impose the instantaneous ap- 
proximation onto the kernel as Vp (p, k) = V(p, k). Tak- 
ing the concrete steps given in Ref. [18j |. we can obtain 
this kernel by a Fourier transformation of the Cornell 
potential which contains a linear piece and a Coulomb- 
type piece, —W{r) — I ® IV s (r) + 7^ ® r y fl V v (r), where 
V s (r) = Xr + Vo and V v (r) = — g^t. To avoid any in- 
frared divergence, a convergent factor e _/3r is introduced 
into the potential, 



V s (r) = ±(l-e-P r ) + V , V v {r) = -\^ 
p 3 r 



-0r 



(8) 



After the Fourier transformation, the kernel in momen- 
tum space reads [3: -iV(q) = I ® IV s (q) + 7^ ® 
7m K(q), where 



V s (q) = -(\/p + V )5 3 (q) + 



A 1 

^2 (q 2 +/3 2 ) 2 



K(q) 



2 a a (q) 
'3tt 2 (q 2 +I3 2 ) 



(9) 



and the effective coupling constant is given by a s (q) 

2 

4tt/9 log (a + — ) with a being a parameter which 



'«CD 



freezes the running coupling constant at low energy. 
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With this kernel, the BS equation (J5]) is written as 
i 1 + i 



x P (p) 



fe. + E /2 + is) (p e -Eo/2- ie) 2 

""" k(qk P (fc) + K(q)7 M X P (fc)7M 



(2tt) 4 



1 



(10) 



where q = p k. Substituting the BS wave functions 
(O and ([7]) into Eq. (fit?]) we get the following component 
equations for vector and pseudo-scalar quarkonia: 



fa(b)ip) = 



(p e +E Q /2 + ie){p t -E /2-ie) 

^(K-K)(q)/ Q(6 )(fc). (11) 



2*) 



For later convenience, here we also write out the BS equa- 
tion for the conjugate wave function, 



XpO) 



(ft + £o/2 + ie)fe -Eb/2-ie) 2 

d 4 fc r 



(2tt) 4 

1 + jl 



V s (q)x p (k) + V v {c L )rX P {k) 1 ^ 

(12) 



To solve the BS equation (flTj) . for convenience, we can 
choose a coordinate frame in which the binding system 
is static. The BS wave functions in this frame are given 
in Eqs. Q ([7]), where f a (p) and fb(p) are Lorentz scalar 
functions. Then the longitudinal component p e = p . 
After carrying out the integrations over p° and k° on 
both sides of the BS equation along some proper contour, 
Eq. (fTTj) becomes 1 



/a(6)(p) 



d 3 k 

E J Wf 



(K-^)(q)/ Q(b) (k), (13) 



where we have defined the instantaneous wave functions 

by / (5)(P) = / ^-fa(b)(p)- 

Normalization of the BS wave functions. 

The normalization condition of the BS wave function \ p 
for a vector meson can be written as 



d 4 pd 4 p' 
(2tt) 8 

- 2P°<W , 



x p s \p) 



d 



Ip a (p,p') 



x { ;'\p') 



(14) 



1 In pratical calculation, to obtain the reasonable BS wave func- 
tions, we will replace Eq by M a ^ — 2^k 2 + m.g in the following 

BS equation. This replacement is equivalent to regaining some 
part of 1/mQ effects in the calculation. 



where I P [p,p') = r (27r) 4 <5 4 (p - p') S~ 1 (pi)S~ 1 (p 2 ) , s 
and s' are the spin indices of the vector meson, and P° = 
\J~P'a + Ma is the energy of the bound state. Multiplying 
by 5 ss i on both sides and summing over s and s' , one has 
the following normalization equation (in the static frame 
of the meson): 



Sv° a f d 3 p ~ ~ 2 



E 



(15) 



'0 J (27T) 

for the vector meson. f a \ and f a 2 are defined as follows 
d 3 k 



/ol(p) 



M„ 



/a2(p) = M a 



(2tt) 3 

rf 3 k 

(2tt) 3 



y s (p-k)/ (k), (16) 



K(p-k)/ a (k). (17) 



From the BS equation (JT3J), /oi - f a 2 = -E M a f a , one 
can see that the normalization equation (|15p can be writ- 
ten as 



4M a 



rf 3 p 

(2tt) 3 



/a( P ) 2 = F 



(18) 



The normalization equation for the pseudo-scalar meson 
has completely the same form as Eq. (fT5|) (with the sub- 
script "a" replaced by "6"). 



III. DOUBLE QUARKONIUM PRODUCTION 
FROM e+e" ANNIHILATION 



Now we turn to discuss the exclusive process in the 
electron-positron collisions with the final state of two 
heavy quarkonia. The relevant Feynman diagrams for 
the process e + e _ — > J/i/jr] c are depicted in Fig. [T] Since 



Pi e 



P2 e 




J/i> Pa (9a) 



Vc Pb (?(,) 



+ permutations 

FIG. 1: Feynman diagrams for double charmonium produc- 
tion from electron-positron annihilation. 

the s-channel gluon turns into cc in this process, hence 
the gluon is hard and the energy scale for this process is 
large. Therefore, the perturbative calculation in QCD to 
lowest order of a s expansion is expected to be sufficient. 
We will not consider the contributions from diagrams of 
higher-order in a s . Furthermore, we do not consider any 
higher order corrections from the pure electromagnetic 
interactions, which should be small and were discussed 
in Ref. @. 
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One of the amplitudes in FigQ]can be written as 

g^ee^ fd*q a d% 
Mi = C - — v{p 2 )^^u{pi) 



xTr 



1 



X P „ 7/3 ■ 



(2tt) 8 

9^ 
k 2 ' 



(19) 



where C = 4/3 is the color factor, y/s is the total energy 
in the center-of-mass frame, and eg is the electric charge 
of the quark Q. The conjugation of the BS wave function 
is defined by x P = 7°(Xp)^7°- The momenta of the quark 
and anti-quark within the final state are q[ = \P a + <Za , 

?3 = \Pa - q a , 94 = \Pb + qb, q2 = \Pb - qb , and 
the momenta in the gluon and the quark propagators are 
given by 



k -- 
3i 



93 + <74 
= q[+k 



\{P a + Pb)-q a + qb, (20) 



Pa + ~ 2 P b 



qb , 



(21) 



respectively. Since an integration is needed to obtain the 
amplitude (TT9"]) and since the propagators of the quark 
and the gluon depend on the relative momenta q a and 
qb one can expect that the calculation is very compli- 
cated. To simplify the calculation, we assume that the 
propagators of the quark and the gluon are independent 
of relative momenta q a and qb (see, e.g. Ref. [21|). This 
simplification is appropriate since the masses of heavy 
quarks are large compared with the relative momenta, 
which are of order a s m,Q. Then the momenta q\ and k 
of the propagators are large compared with the relative 
momenta q a and q b . One may expect that, in the heavy 
quark limit, the calculation without taking into account 
the relative momenta should be exact. 2 

With the above approximation in mind, the momenta k 
and gi can be written as: k » (P a + Pb)/2, q\ P a +Pb/2 
which lead to k 2 » s/4, q\ » s/2 + A/ 2 /2-M 2 /4 » s/2 + 
m 2 . Furthermore, we will make use of the approximation 
Mji^ w M Vc w 2m c in the calculation throughout this 
paper. Then the amplitude (|19|) can be written as 



Mi 



g s e v{p2)-l^u{ Pl ) 



d 4 q a d 4 q b 



3 2 s 3*s w ' H 1 sru J ( 2?r )4 ( 27r )4 
xTr [x^Hqaha (4i + m c ) 7 ^X p >) 7 Q ] ■ (22) 



From the component expressions of the BS wave func- 
tions © and 0, we can see x P ^ = ^ (s) (l + ^ )M / and 
X P = 75(1 + tfb)Mbfb' 3 Substituting the BS wave func- 



2 The energy scale fj, is of the same order of ttiq and then, in the 
heavy quark limit, d s (/t) ~ 1/ l°g( m Q/^QCD) - * 0. 

3 We have rotated the phase of the wave function to make f a and 
ft be real. 



tions into the trace in the above amplitude, one achieves 
2 14 7r 2 a s a Rm 



Mi = -i- 



T7l c C^vpo- 



(s)u pp per 
° 1 a 1 b 



x v(p 2 h fi u(pi)il) a ip b . 



(23) 



where a s = g 2 /4ir, a cm = e 2 /47r, ip a and ip b are two 
numbers defined by the integrations over f a and /&, re- 
spectively, 



d 4 q 



d 3 q 



<Pa(b) = J ^4 /«(6)(?) - J 7^3 /a(6)(q) • (24) 

The total amplitude for the process e + e~ — > J/iprj c can 
be obtained by summing over all the amplitudes of the 
diagrams shown in Fig. [1] 

The unpolarized total cross section (see e.g. Ref. |22 |) 
is obtained by summing over various J / ip spin-states and 
averaging over those of the initial state e + e~, 



1 yfs — 16m 2 



32tt 



,3/2 



/i^lMataifdcosfl, (25) 



where the masses of the electron and positron are ignored 
in the calculation. The explicit expression for the total 
amplitude |.Mtotai| 2 , which is the sum of all diagrams 
shown in Fig. [TJ is written as 



total 



- — m c {—62m c +t +u ) 



spins 



3 4 



(26) 



where t = (pi — P a ) 2 and u = (pi — Pb) 2 are the Mandel- 
stam's variables. 

Numerical results. 

The parameters in the calculation will be taken to be 
0Gj|: a = 2.7183, = 0.06 GeV, a om « 1/137, a s = 
0.26, A= 0.2 GeV 2 , A QCD = 0.26 GeV, m c = 1.7753 
GeV. In the interaction kernel Vo = 0.415 GeV for J/tp 
and Vo = 0.525 GeV for r/ c . With these parameters, 
we can solve the BS equations numerically and the wave 
functions f a and f b are plotted in Fig(2J 
Then we have 

tp a = 0.1020 GeV , Vfc = 0.1037 GeV. (27) 

Consequently, the total cross section is obtained as 

cr(e+e- -> J/ip rjc) = 22.3 fb (28) 

If we vary m c and a s within 10% we can get an errors 
±1.6/6. 

The above analysis for the exclusive process e + e~ —> 
J I tp r\ c can be applied, with only a little modifications, to 
the exclusive process e + e~ — > T(lS)rjb- ?7fc(9434) is the 
lowest-lying pseudo-scalar bb state (for discussions about 
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CO 0.5 

CO 



FIG. 2: BS wave functions for J/tp and 7/ c in the heavy quark 
limit, /„ (solid line) and /t, (dashed line), respectively. 



i]b (9434), see e.g. Ref. [23j). The mass of the b quark is 
m b = 5.224 GeV [H[. In the interaction kernel V = 0.62 
GeV for T(15) and V = 0.64 GeV for rj b . With these 
parameters, we can solve numerically the wave functions 
f a and f , which lead to 

ip a = 0.1123 GeV , ip b = 0.1124 GeV. (29) 

Then the cross section is predicted to be 

a(e+e- -> T{1S) r) b ) = (0.16-0.06) fb (30) 

for the range of the total energy y/s — (25-30) GeV. 

IV. CONCLUSIONS AND DISCUSSIONS 

In this paper, we study the exclusive processes of e + e~ 
annihilating into two quarkonia in terms of the BS for- 
malism. We find that, in the heavy quark limit, the 
cross section is a[e + e~ — > J/tprjc] = 22.3 fb, which is 
compatible with the Babar's data, cr[e + e~ — > J/'ipr/c] > 
17.6 ± 2.8 ± 2.1 fb 2], and the Belle's data, cr[e+e~ -> 
J/tpVc] > 25.6 ± 2.8 ± 3.4 fb Q. Because the BS for- 
malism is established based on the relativistic quantum 
field theory, one has a strong reason to believe that some 
(perhaps not all) relativistic effects are automatically in- 
cluded in the calculations. The missing part may come 
from the instantaneous approximation which is necessary 



for solving the BS equation. Since the approximation is 
proved to be reasonable in theoretical calculations for 
other similar processes, we may be convinced that the 
missing part is not significant. Thus we expect that the 
non-leading-order contributions, from extra 1/rriQ which 
is indeed a relativistic effect, and a s corrections, should 
be small. Our result is different from those given in 
Refs. @, 0, El O in the NRQCD framework, where the 
leading-order contribution is too small to be compara- 
ble with the experimental data. In order to reduce the 
discrepancy between the experimental data and the the- 
oretical predictions based on NRQCD, large non-leading- 
order contributions (including a s or/ and relativistic cor- 
rection^)) are required in the NRQCD framework and 
the value of the total non-leading-order correction is 
nearly twice of the leading-order one. 

We also calculate the cross section for the exclusive 
process with two bottomonia as the final state, <r[e + e~ — > 
T(lS)r)b] — (0.16-0.06) fb corresponding to the range of 
the total energy y/s = (25-30) GeV. Since m b ^> m c , 
the non-leading-order contributions, including effects of 
higher orders in l/m b and a s (2m b ) expansions, should be 
much smaller than those in the charmonium case. There- 
fore, we expect that the calculation in the bottomonium 
case is much more precise than that in the charmonium 
case. Even though such processes cannot be measured 
at present due to the constraint of the available energy 
range at the B-factories, the planned ILC will be a pow- 
erful facility to testify this result. 

Due to the obvious advantage of the BS formalism for 
dealing with the processes where heavy flavors are in- 
volved, we may hope that the obtained results are close 
to reality. 

It will be interesting to study the non-leading-order 
contributions in the formalism used in this paper and 
check whether our expectation about the non-leading- 
order contributions is right. This task is beyond the scope 
of this paper and will be discussed elsewhere. 
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